Anomalous heat conduction and anomalous diffusion in one dimensional systems 
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We establish a connection between anomalous heat conduction and anomalous diffusion in one 
dimensional systems. It is shown that if the mean square of the displacement of the particle is 
(Aa; 2 ) = 2Dt a (0 < a < 2), then the thermal conductivity can be expressed in terms of the system 
size L as k = cL 13 with (3 = 2 — 2/ a. This result predicts that a normal diffusion (a = 1) implies 
a normal heat conduction obeying the Fourier law (/3 = 0), a superdiffusion (a > 1) implies an 
anomalous heat conduction with a divergent thermal conductivity (/3 > 0), and more interestingly, a 
subdiffusion (a < 1) implies an anomalous heat conduction with a convergent thermal conductivity 
(f3 < 0), consequently, the system is a thermal insulator in the thermodynamic limit. Existing 
numerical data support our results. 



Does heat conduction in one dimensional (ID) sys- 
tems obey the Fourier law? If it does, what are the 
necessary and sufficient conditions? If it does not, then 
what is the reason and how the thermal conductivity di- 
verges/converges with the system size LI These ques- 
tions have attracted increasing attention in recent years 
[1-21]. Although some progress have been achieved, 
many puzzles remain. For example, in an attempt to 
establish a connection between heat conduction and the 
underlying microscopic dynamics, there exist some con- 
troversial examples. In the ding-a-ling model Casati et al 
[2] show that as onset of global chaos the heat conduction 
crosses over from an abnormal one to a normal one obey- 
ing the Fourier law. It is thus concluded that the chaos is 
a deciding factor. Later on, in order to show the exponen- 
tial instability is a necessary condition, Alonso et al [12] 
studied the heat conduction in a Lorentz gas channel, a 
quasi ID billiard with circular scatterers, and found that 
the heat conduction obeys the Fourier law. However, the 
results from ID Ehrenfest gas channels [15], in which the 
Lyapunov exponent is zero, show that the Fourier heat 
law might not have any direct connection to the under- 
lying dynamical chaos, because the heat conduction can 
be normal and abnormal, depending on whether or not 
the disorder is introduced. 

Recently, a quasi ID triangle billiard model, which con- 
sists of two parallel lines of length L at distance d and 
a series of triangular scatterers, has been introduced and 
studied [16]. In this model, no particle can move between 
the two reservoirs without suffering clastic collisions with 
the triangles. Therefore this model is analogous to the 
Lorentz gas channel studied in [12] with triangles instead 
of discs and the essential difference is that in the trian- 
gular model the dynamical instability is linear and there- 
fore the Lyapunov exponent is zero. It is found that the 
motion inside the irrational triangle channel (the inter- 
nal angles are irrational multiples of n) is diffusive and 
has a normal heat conduction. Therefore deterministic 



diffusion and normal heat transport which are usually as- 
sociated to full hyperbolicity, can take place in systems 
without exponential instability. Another example is the 
Fermi-Pasta-Ulam (FPU) model [4] which has non-zero 
Lyapunov exponent, however the heat conduction in this 
model does not obey the Fourier law. 

The heat conduction in the rational triangle model (the 
internal angles are rational multiples of 7r) and in the 
FPU model is anomalous and does not obey the Fourier 
law, the thermal conductivity k diverges with system size 
L as L@ with = 0.22 for the rational triangle model 
[16], and 0.34 < < 0.44 for the FPU model [4]. Indeed, 
similar divergent behavior has been observed in many 
ID systems. For example, in the binary hard sphere 
model [19,20], 0.22 < (3 < 0.35, in single wall nanotubes 
0.22 < < 0.37 [18], and in many classical lattices such 
as the harmonic lattice, (3 = 1 [22], disordered harmonic 
lattice, 0=1/2 [23], and the Frenkel-Kontorova (FK) 
model under the condition of T/K » 1, < < 1 [14], 
where T is temperature, and K is the effective amplitude 
of a sinusoidal on-site potential. 

Obviously, a universal value of (3 does not exist, it 
differs from model to model. Most recently, Narayan 
and Ramaswamy [21] show theoretically that in a ID 
momentum-conserving continuous system, the heat con- 
duction is anomalous, and the thermal conductivity di- 
verges with system size L as L 1 / 3 . Up to now, in all avail- 
able numerical results only the heat conduction in a (5,5) 
single wall nanotubes [18] shows an exponent {(3 « 0.32) 
close to this 1/3 [24]. Despite the conduction mechanism 
is similar, a (10,10) single wall nanotube shows different 
value [18] for unknown reasons. The numerical results 
from other models such as the FPU model, the harmonic 
model and other billiards models deviate largely from 
this value for reasons to be investigated. 

On the other hand, even if the momentum conservation 
breaks down, the heat conduction can be anomalous such 
as that in the Frenkel-Kontorova model [14]. The ques- 
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tion becomes how to explain this anomalous heat con- 
duction, in particular the value of the exponent in the 
thermal conductivity. A general theory is still lacking. 
The only existing theory is for the ID harmonic chain 
[22], in which the phonons transport along the chain bal- 
listically and the thermal conductivity, k, diverges as L, 
i.e. (3=1. 

In this Letter, we would like to find a microscopic ori- 
gin of the anomalous heat conduction observed in many 
ID models. We shall not restrict to any specific model. 
This should give us a more general way to understand 
the heat conduction in ID systems. 

As is well known that, depending on the value of ex- 
ponent a in the mean square of displacement of the par- 
ticle, (Aa; 2 ) = 2Dt a with < a < 2, ID microscopic 
motion can be classified into ballistic motion, a = 2, su- 
pcrdiffusion, 1 < a < 2, normal diffusion, a = 1, and 
subdiffusion a < 1. Ballistic transport is observed in 
the harmonic lattice. Normal diffusion shows up in the 
FK model in a certain parameter regime [5], the disor- 
dered FPU model [10], the Lorentz gas channel [12], the 
disordered Ehrenfest gas channel [15], the irrational tri- 
angle channel [16], and the alternative mass hard-core 
potential model [17]. In some billiard models, superdif- 
fusion [15,16,25-27] and subdiffusion [25] are observed. 
Supcrdiffusion and subdiffusion can be studied from the 
fractional Fokker-Planck equation, for detailed theoret- 
ical investigation and discussion about the anomalous 
diffusion, please refer to review articles [27,28] and the 
references therein. 

To establish a connection between the microscopic pro- 
cess and the macroscopic heat conduction, let's consider 
a ID model of length L whose two ends are put into 
contact with thermal baths of temperature Tl and T R 
for the left end and the right end, respectively. Sup- 
pose the energy is transported by energy carriers (they 
are phonons in lattices and particles in billiard channels) 
from left heat bath to the right heat bath and vice versa. 
If the mean square of displacement of the carrier, with 
velocity v, inside the system can be described by, 



(Aa; 2 ) = 2Dv a t a 



(1) 



then the so-called "mean first passage time" (MFPT) is 
[29] 



(t LR ) = 



(2) 



Obviously, if the Id system is isotropic, the MFPT for 
the carrier travelling from the right to the left end (t R L) 
is the same as (tLR)- 

If the heat bath is a stochastic kernel of Gaussian 
type, namely, the probability distribution of velocities 
is p{v,T) = 47rw 2 exp(-w 2 /2T)/(27rT) 3 / 2 , the MFPT be- 
comes, 



{tLR) = 
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(3) 



We define the heat current as the energy exchange be- 
tween two heat baths in unit time. Thus the current 
induced by a carrier (m = 1) with velocity v moves from 
left to right and comes back is: 



._ J °°i (p(v, Tl) - p{v, T R )) = T L - T R 

3 " (t L R) + (tRL) 2(t LR ) ' 



(4) 



If the temperature difference between the two baths is 
sufficient small so that VT = (T R - T L )/L, then the 
thermal conductivity, k = —Lj/VT, is 



n = cL f) , = 2-2/a, 



(5) 



and the constant, c = 3nV2^a{irV~D/2) 2 / a VT/(32~f). 

Eq. (5) is the central result of the paper [30]. It con- 
nects heat conduction and diffusion quantitatively. The 
main conclusion is that an anomalous diffusion indicates 
an anomalous heat conduction with a divergent (conver- 
gent) thermal conductivity. More precisely, our result 
tells us that: a ballistic motion means thermal conductiv- 
ity proportional to the system size L, a normal diffusion 
means a normal heat conduction obeying the Fourier law, 
a supcrdiffusion means a divergent thermal conductivity, 
a subdiffusion means a zero thermal conductivity in the 
thermodynamic limit. In the following, we compare our 
results with the existing analytical and numerical results. 

A ballistic motion, a = 2, leads to a divergent ther- 
mal conductivity k cx L. The only existing analytical 
result is heat conduction in a ID harmonic lattice. It 
is known that heat is transported by phonons in lattice 
model. Because there are no resistance and umklapp 
process, the phonons transport ballistically in harmonic 
lattice model, thus a = 2. From our formula (5), the 
thermal conductivity in the ID harmonic lattice diverges 
as LP with = 1, this is exactly what was shown by 
Lebowitz et al [22] ("*" in Fig. 1). 

A normal diffusion, a = 1, means that the thermal 
conductivity is a size independent constant, = 0, i.e. 
the heat conduction obeys the Fourier law. For example 
in the ID Frenkel-Kontorova model [5], in a certain range 
of parameter such as T/i-T << 1 [31], the phonons trans- 
port diffusively [7] , thus the thermal conductivity is finite 
and independent of the system size L. The disordered 
FPU model also has a finite thermal conductivity due 
to the random walk like scattering process in the chain 
[10]. Other ID models showing normal diffusion and nor- 
mal thermal conduction are: the ID Lorentz gas channel 
[12], the ID disordered Ehrenfest gas channel [15], the 
ID irrational triangle channel [16], the alternative mass 
hard-core potential model [17], and some ID polygonal 
billiard channels with certain rational triangles [25] . "*" 
in Fig. 1 represents all models with normal diffusion. 
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FIG. 1. The a — [3 plot. Normal diffusion: * represents 
models with a normal diffusion and a normal heat conduction, 
i.e. a = 1 and (3 = 0, such as the Lorentz gas channel[12], the 
Frenkel-Kontorova model [5], the <j> 4 model [9], the disordered 
FPU model[10], the disordered Ehrenfest gas channel[15], the 
irrational triangle channel [16], the alternative mass hard-core 
potential model[17], and some rational polygonal channel[25] 
etc. Ballistic motion: * represents the ballistic transport, i.e. 
a = 2 and (3=1, such as the ID harmonic lattice model. 
Super-diffusion: y, ID Ehrenfest gas channel with right an- 
gle triangle scatterers[15]; O, ID channel with rational tri- 
angle scatterers[16]; A, polygonal billiard channel with one 
irrational ((v5 — l)7r/4) and one rational (tt/3) triangle; O, 
a ID triangle-square channel[26]. Subdiffusion: the polygo- 
nal billiard channel with one irrational angle ((a/5 — l)7r/4) 
and one rational angle (7r/4)[25], □, from the channel length 
1 < L < 40; >, from the channel of length 40 < L < 80. 

A superdiffusion, 1 < a < 2, implies an anomalous 
heat conduction with a divergent thermal conductivity 
L 13 . The exponent < (3 < 1 differs from model to 
model. Here we take the billiard models as our examples 
because they are very clean, and both the diffusion and 
thermal conductivity in these models can be calculated 
very accurately. The first example is the ID Ehrenfest 
gas channel in which the scattering obstacles are isosceles 
right triangles periodically post along the channel [15]. In 
this model one has a — 1.672. From our analytical result 
(5), the thermal conductivity should diverges as L@ with 
(3 = 2 — 2/q = 0.804 which agrees with the result from 
simulation of heat conduction (3 — 0.814 [15] ("v" m Fig. 
1). The second example is the ID channel with triangles 
whose inter angles are rational multiples of tt [16]. This 
model shows a superdifussion with a = 1.178. The di- 
vergent exponent of thermal conductivity is (3 = 0.302. 
This exponent is slightly larger than the one obtained 
from thermal conductivity simulation (3 — 0.22 ("O" m 
Fig. 1). This deviation is due to the finite size effect in 
the heat conduction simulation. 

A subdiffusion, a < 1, results in an anomalous heat 
conduction with a convergent thermal conductivity, i.e. 
k oc L 13 , with (3 < 0. This is an interesting result im- 
plying that the system becomes a thermal insulator in 
the thermodynamic limit L — ► oo. Although there are 



many examples showing subdiffusion [32-36], a system- 
atic study on the heat conduction in such kind of systems 
is still lacking. The only existing example is the heat 
transport in a polygonal billiard which supports our con- 
clusion [25] . Most recently, Alonso et al [25] show that in 
a very special configuration, a = 0.86, and the thermal 
conductivity goes as k ~ £-°- 63 ("□" m Fig. 1). As L 
goes to infinity, the thermal conductivity goes to zero. 
According to our formula (5), if a = 0.86, (3 = —0.33 
which is larger than the one obtained by Alonso et al [25] . 
This is not a surprise, because the channel length in their 
study of thermal conductivity is too small (L < 40). If 
the channel is longer, the value of (3 will become much 
more closer to our theoretical estimation ((3 = —0.33). 
To demonstrate this, we extend the thermal conductiv- 
ity simulation from L E [1,40] used by Alonso et al [25] 
to L £ [40,80], and find that (3 = -0.48 (">" in Fig. 1) 
which is more closer to (3 = —0.33 than the one obtained 
by Alonso et al. If L — ► oo, one can expect (3 goes to 
-0.33. 

All numerical results are summarized and represented 
in Fig 1, where we draw (3 versus a, and compared with 
Eq. (5). As is shown that, Eq. (5) is exact for both nor- 
mal diffusion and the ballistic motion. The agreement 
with most existing numerical data is good. However, 
discrepancies remain for some models mainly due to the 
limited numerical simulations. The best data close to 
curve (3 — 2 ~ 2/ais the simulation from ID Ehrenfest 
gas channel [15]. This is because the channel length used 
in the simulation is the longest one (L ~ 10 3 ) among all 
the available data. 

In summary, we have established a connection between 
anomalous heat conduction and anomalous diffusion in 
ID systems. Our central result Eq. (5) includes all possi- 
ble cases observed in different classes of ID models, rang- 
ing from subdiffsion, normal diffusion, and superdiffusion 
to ballistic transport. Several conclusions can be drawn: 
(1) A normal diffusion leads to a normal heat conduction 
obeying the Fourier law. (2) A ballistic transport leads 
to an anomalous heat conduction with a divergent ther- 
mal conductivity k cx L. (3) A superdiffusion leads to 
an anomalous heat conduction with a divergent thermal 
conductivity in thermodynamic limit. (4) More impor- 
tantly, our result predicts that a subdiffusion system will 
be a thermal insulator. Existing numerical data support 
our results. 

We should mention that the subdiffusion process has 
been observed in many real physical systems such as 
highly ramified media in porous systems [32], percola- 
tion clusters [33], exact fractals [34], the motion of a 
bead in a polymer network [35] , charge carrier transport 
in amorphous semiconductors [36]. Any numerical simu- 
lation and/or real experimental measurement of thermal 
conductivity in these systems will be very interesting and 
will allow one to test the theory given in this Letter. More 
importantly, it will have a wide application in designing 
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novel thermal devices. 
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